Nondiffusive suprathermal ion transport in simple magnetized toroidal plasmas 



O 

(N 

o 

Q 

oo 



K. Gustafson, P. Ricci, I. Furno, and A. Fasoli 

Ecole Polytechnique Federate de Lausanne (EPFL), Centre de Recherches en Physique des Plasmas, 
Association Euratom- Confederation Suisse, CH-1015 Lausanne, Switzerland 
(Dated: December 9, 2011) 

We investigate suprathermal ion dynamics in simple magnetized toroidal plasmas in the pres- 
ence of electrostatic turbulence driven by the ideal interchange instability. Turbulent fields from 
fluid simulations are used in the non-relativistic equation of ion motion to compute suprathermal 
tracer ion trajectories. Suprathermal ion dispersion starts with a brief ballistic phase, during which 
particles do not interact with the plasma, followed by a turbulence interaction phase. In this one 
simple system, we observe the entire spectrum of suprathermal ion dynamics, from subdiffusion to 
superdiffusion, depending on beam energy and turbulence amplitude. We estimate the duration 
of the ballistic phase and identify basic mechanisms during the interaction phase that determine 
the character of suprathermal ion dispersion upon the beam energy and turbulence fluctuation 
amplitude. 

PACS numbers: 52.65. Cc, 52.25.Fi, 52.30.Ex, 52.35. Ra, 52.40.Mj 
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Most plasmas are characterized by the presence of 
suprathermal particles, possibly generated by turbulent 
acceleration, external sources, or, as in the case of fu- 
sion devices, nuclear reactions. Understanding the ba- 
sic phenomena that determine the suprathermal particle 
dynamics is a key challenge for the description of a wide 
range of plasma systems, ranging from magnetically con- 
fined plasmas for fusion [iH^I to space plasmas [Hj. 

In this Letter, we theoretically characterize suprather- 
mal ion dynamics in the simple magnetized torus (SMT) 
configuration Q . An SMT confines a plasma with a verti- 
cal magnetic field, Bv, superimposed on a toroidal mag- 
netic field, B^, creating helicoidal field lines that ter- 
minate on the vessel. This configuration is of interest 
to the plasma turbulence and fusion communities since 
it offers a simple and well-diagnosed testbed in which 
to study the basic physics of plasma turbulence and the 
associated transport of heat and particles, allowing pa- 
rameter scans that are not possible in more complicated 
configurations. In the SMT, turbulence has been charac- 
terized through global simulations validated against 
experimental data [1, These are unique simulations, 
which evolve the plasma dynamics resulting from the in- 
terplay between the plasma source, losses at the vessel, 
and turbulence, with no separation between equilibrium 
and fluctuating quantities. They provide the turbulent 
fields for integrating suprathermal ion trajectories. 

The SMT incorporates, in a simplified form, the fun- 
damental elements determining suprathermal ion dynam- 
ics, specifically: the Larmor gyration, the drifts related 
to the curvature and gradient of the magnetic field, the 
E X B drift, and the polarization drift. The latter two 
are related to plasma turbulence and strongly dependent 
upon its topological properties. The relative simplicity 
of the SMT allows comprehensive quantification of the 
interplay between these phenomena. Since the key ele- 
ments are the same, the framework established here can 
be applied for interpreting suprathermal ion dynamics 



in more complicated and diverse contexts. Examples in- 
clude fusion devices with high energy neutral beams and 
a-particle production, cosmic ray propagation, and solar 
wind interaction with the magnetosphere. 

The primary diagnostic for studying the dispersion of 
suprathermal ions is the variance <j\{t) = {^R^) ^ t"''^ of 
their radial displacements, 6R = R{t) — i?(0), where <> 
is an ensemble average over many particle trajectories. 
By numerically integrating the trajectories of suprather- 
mal ions in simulated SMT turbulent fields, and by ex- 
ploring wide ranges of particle energy and turbulence am- 
plitude, we show that the ions have a complex motion, 
which in general cannot be considered diffusive. Our 
simulations show that suprathermal ion dispersion starts 
with a brief ballistic phase, during which particles do not 
interact with the plasma, resulting in 7^ ~ 2. This phase 
is followed by a turbulence interaction phase, which sur- 
prisingly shows the entire spectrum of suprathermal ion 
spreading: superdiffusive {jr > 1), diffusive (7^ = 1), 
or subdiffusive (7^,; < 1), depending on particle energy 
and turbulence amplitude. We provide an estimate of 
the duration of the ballistic phase and we identify the 
mechanisms that determine, in the interaction phase, the 
dependence of "fn upon the beam energy and turbulence 
fluctuation amplitude. 

For an SMT exam ple, we refer to the parameters of 
the TORPEX device [i3,|iH, namely B„ < B^, /3 < 1 
and Ti -C . In TORPEX, a localized source of plasma 
on the high-flcld side of the torus is generated by mi- 
crowave absorption at the electron-cyclotron and upper- 
hybrid resonances. A suprathermal ion beam is provided 
by a miniaturized lithium ion source [l^ . A num- 
ber of turbulence regimes have been characterized both 
experimentally and theoretically 0] for TORPEX. Here, 
we focus on the ideal interchange instability, which is 
dominant for a sufficiently high By and low plasma re- 
sistivity. In this regime, fc|| = 0, and the wavelength of 
the dominant mode along z (the direction perpendicu- 
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FIG. 1. (Color online) Suprathermal ions in a schematic of 
TORPEX. A section of the torus vessel is indicated, as well as 
a magnetic field line (dashed black line) . Three suprathermal 
ion trajectories (red, green and blue solid lines) for a, £ — 20 
injection are shown. The injection cone is indicated in red, 
while the black arrows indicate several ions at tQ — 130. A 
snapshot of the electrostatic potential, $, obtained from the 
simulations in Ref. 01 and used to integrate the ion trajecto- 
ries, is shown on two poloidal cross-sections. 



lar to both B and the radial direction, which is tilted 
by an angle 9 — tan^^ {Bi,/B^) <^ 1 with respect to 
the vertical direction), is given by the return length of 
the field line on the poloidal plane, i.e. \z = A, where 
A = 2TTRoBy/ Bel), with Rq being the SMT major radius. 
Since fcy = 0, turbulence in this regime can be described 
by two-dimensional simulations in the tilted plane per- 
pendicular to the magnetic field. 

We integrate the non-relativistic ion equation of mo- 
tion with the Lorentz force in the SMT magnetic con- 
figuration, considering the time-dependent electric field 
provided by the simulation described in Ref. 0]. An ex- 
ample of the simulated electrostatic potential for TOR- 
PEX, $, is shown in Fig. [T]in two perpendicular planes 
at different toroidal locations. Two distinct poloidal re- 
gions with different features of plasma turbulence have 
been identified previously in both simulations and exper- 
iments. First, a mode region exists on the high-field side 
of the SMT, where the plasma is generated and a coher- 
ent ideal interchange mode is present. Second, a region 
on the low-field side |13j with lower plasma density and 
temperature is marked by intermittent structures termed 
blobs [li,[l5|. 

Injection conditions for the ions are chosen to mimic 
the suprathermal ion source in TORPEX, focusing on 
parallel injection, such that the axis of the injection cone 



is directed along a field line using a Gaussian angular 
distribution with a O.lrad variance. The initial velocities 
have a Gaussian distribution with a spread cr^g = O.lwg, 
with vq being the mean initial velocity, injected at the 
point R = Rq, a region where turbulence is transitioning 
from the mode region to the blob region. We consider 
the ions as tracer particles, such that they do not in- 
fiuence background fields, and ignore suprathermal ion 
collisions. These two assumptions are motivated by the 
experimental conditions in TORPEX. 

We now detail the main elements determining a 
suprathermal ion trajectory in SMT turbulent fields. In 
the parallel direction, particle velocity is essentially unaf- 
fected since fcy = 0. In the perpendicular plane, ion tra- 
jectories result approximately from the combined effects 
of four elements. These are the gyromotion, with Lar- 
mor radius p = v±/Q and Larmor frequency f2 = qB/m 
(in and q are the suprathermal ion mass and charge, re- 
spectively), the drifts related to the curvature and radial 
gradient of the magnetic field, wvb, the E x B drift, and 
the polarization drift. 

The drift velocity uvb = (^^1/2 + v'pez/(flR) is ori- 
ented purely in the z direction and dominated by the 
curvature term for parallel injection. For the E x B 
drift, we note that the time- averaged electric field has a 
radial component that causes a drift in the z direction. 
The fiuctuating electric field leads to alternating displace- 
ments in the radial and z directions. The size of these 
displacements is determined by the size and amplitude 
of the fluctuating vortex and blob-like structures in the 
turbulence, and by the Larmor radius. In fact, the E x B 
drift accounts for the gyroaveraged electric fleld. In the 
case where the suprathermal ion Larmor radius is signif- 
icant compared to the scale of the turbulence, ^aP ^ 1 
(/cA = 27r/A), the gyroaverage decreases the magnitude 
of the E X B velocity with respect to cases for which 
f^AP ^ 1- Finally, the polarization drift causes modest 
but observable particle energization, similar to that seen 
in Ref. [l^. This can affect the transport slightly by 
increasing the average Larmor radius. 

In the SMT, the turbulent E x B drift is thus the 
sole cause of suprathermal ion beam radial dispersion, 
which we quantify by the variance cr|j(f) of the radial 
displacements. Figure [5] illustrates an example of the 
time evolution of o-|j(i) with the three phases that appear 
in the simulations. These phases are characterized by 
different values of the 7ij exponent, as determined by 
fitting a line in log-log plots of (t) . The phases can be 
categorized as ballistic, interaction and asymmetric, as 
shown in Fig.[5J They are distinct, but tend to transition 
smoothly from one into the next, as measured by the 
dispersion exponents. 

The early ballistic phase is the relatively brief period 
with '-^ 2 before the ions interact significantly with 
the turbulence and magnetic field. In fact, the sim- 
ulations show that the ballistic phase lasts until t ~ 



iLnm{2/V,,mv±/{qE)}. In principle, a gyrocenter bal- 
listic phase may be present, during which the particle 
dispersion shows 7^ ~ 2 superimposed on a Larmor os- 
cillation. This gyrocenter ballistic phase would end when 
the gyrocenter speed changes significantly with respect to 
its initial value. In the present context, we estimate this 
phase to last less than one gyroperiod. The simulations 
agree that gyrocenter motion is never ballistic. 

Following the ballistic phase, suprathermal ions inter- 
act with the plasma turbulence and the beam dispersion 
displays nearly constant ^u- This quasi-steady state in- 
teraction phase ends when the ions have spread radially 
enough to sample regions where the local turbulent prop- 
erties are significantly diff^erent than the ones at the injec- 
tion position. At this time, the asymmetric phase begins 
because the beam has spread to sample regions where 
the turbulence properties are significantly different with 
respect to the injection position. Thus, the duration of 
the interaction phase is determined by the spread in the 
radial E x B velocity compared to the width of a roughly 
uniform region of turbulence. 

In Fig. [31 the values of in the interaction phase, as 
obtained from a large number of simulations, are shown 
as a function of the beam energy, which we express in 
terms of £ = S = mwQ/(2Te), and of the amplitude of 
the fiuctuations of the electric field, ^ = e$/Te. Here, 
Tg denotes the electron temperature, averaged over time 
and the z coordinate, and 4> is the root mean square 
fiuctuation amplitude of Both are evaluated at the 
injection point. The fluctuation amplitude is scanned by 
rescaling the fluctuating part of the simulated electric 
field. Due to variations in the plasma properties, this 
is also equivalent to varying ^ and £ by injecting ions 
at different radial positions. At R = Rq, the typical 
experimental value is ^ = 1.5, while the simulations show 
smaller amplitude fluctuations, ^ = 0.8. Figure [3] reveals 
the existence of different regimes for the ion dynamics. 
The complex dependence of 7^ on these parameters is 
now explained. 

Superdiffusive dispersion, i.e. 1 < 7^, is observed in 
a region of the ^ — £ plane, corresponding to k±p 0, 
WVB 0, and ^ > 0.3. In this limit, turbulent struc- 
tures are relatively static with respect to the ions, thus 
allowing ions to move large distances in a single direc- 
tion [l7[. However, if the fluctuations are reduced below 
a certain level, drops dramatically because the ampli- 
tude of the vortex structures is too small for the struc- 
tures to form connections between the center and edge 
of the plasma. This is a topological constraint set by the 
amplitude of the turbulent fluctuations. In Fig. [3J the 
^ ~ 0.3 boundary, for which connected velocity stream- 
lines do not form, is indicated by the horizontal line: jr 
decreases very sharply for ^ < 0.3. 

The boundary marked by the solid curve in Fig. [3] de- 
note the points for which k/^p — 2. Outside this bound- 
ary gyroaveraging gradually decreases the number and 
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FIG. 2. (Color online) Variance of displacements in the ra- 
dial (solid curve), vertical (dashed curve) and parallel (dash- 
dotted curve) directions for £ = 50 and ^ = 0.8 are shown. 
Dispersion exponents 7 are fitted with solid line segments, 
which have error of ±0.1. For the radial dispersion, an initial 
ballistic phase occurs (red-shaded region) with 7i{ ~ 2. This 
is followed by the turbulence interaction phase when 7_r re- 
mains nearly constant. Later, the asymmetric phase (shaded 
blue region) shows an increased value of yR. For the paral- 
lel direction, since there are no forces, 7|| ~ 2 always. The 
z-directed spreading also shows three phases in which the su- 
perdiflFusion is due to i;vb. 

the amplitude of the suprathermal ion radial displace- 
ments. This leads to a reduction of 7^? from superdiffu- 
sive to a diffusive value, 7^ ~ 1. 

The drift velocity wvb has an important effect on the 
suprathermal ion dynamics. Provided that the motion 
along z is sufficiently fast, an effective drift-average of 
the electric field fluctuations reduces the radial disper- 
sion, which is subdiffusive for a significant amount of time 
in the interaction phase. Subdiffusive radial spreading 
occurs if the time required for an ion to traverse a turbu- 
lent vortex along z is significantly smaller than the time 
to traverse the same vortex radially. This time can be es- 
timated as follows. Let tr be the time required to move 
radially across the structure, such that 

tr ^ Lr/vex-b,^, (1) 

where Lr is the radial extent of the vortex, which has 
been estimated [Tsj as Lr ~ Lp/k^. Similarly, we 
define 

Tz ^ LJ {vy-B + VexB.z) (2) 

as the time required to cross a vortex of size Lz ~ 1/fcA 
due to the velocity in z. The curve 

Lr{v^B + VexB.z) 

r — - X (3) 

^zVExB,r 
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FIG. 3. (Color online) Dispersion exponents 7h (colored dots) 
in the interaction phase are displayed in the — £ plane. Error 
on the value of 7_r is ±0.1. Regions for 7_r are demarcated by 
the gyroaveraging condition (A;ap = 2, solid black curve), by 
effective drift-averaging condition (Eq. ((3|, dashed red curve), 
and by the disconnected streamline condition = 0.3, hori- 
zontal line). 



approximately identifies the region where the ions arc 
more likely to complete radial steps before drift-averaging 
makes radial steps less likely. The numerical parameter 
X ~ 5 is empirically observed to correspond to the subdif- 
fusive transition for all values of ^ tested. The condition 
for z;vB averaging, given by Eq. ([3|), is also displayed 
in Fig. [3l confirming the reduction of the suprathermal 
ion dispersion rate to subdiffusive values due to drift- 
averaging. 

Additional simulations with Rq oo confirm the role 
of wvB drift in the suprathermal ion dispersion. In 
these simulations, dispersion is reduced from turbulent 
superdifFusion only by gyroaveraging. The values of 
decrease from 7^ > 1 to 7^ ~ 1 when the Larmor ra- 
dius becomes sufficiently large, but the dispersion never 
becomes subdiffusive. 

Remarkably, we find in one simple system that 
suprathermal ion spreading can be subdiffusive, diffu- 
sive, or superdiffusive depending on the ion energy and 
turbulence amplitude. In previous works, superdifFusion 
and subdiffusion, separately, have been used to model 
plasmas, see e. g. |19l422l |. The coexistence of the three 
regimes in an ad hoc Hamiltonian model was observed in 
Ref. [2^. Nevertheless, a diffusive approach continues to 
be assumed typically. Diffusion may sometimes describe 
suprathermal ion transport over short time and spatial 
scales, since a nondiffusive process may be linearized as 
an "effective diffusion." Our simulations show that this 
approximation is valid only locally, since 'Jr can be dras- 
tically different than unity, and time-dependent. There- 
fore, the effective local suprathermal ion diffusivities can 



show a strong time dependence. In the cases explored 
here, they can be two orders of magnitude away from 
measurements of thermal particle diffusivity computed 
for the ideal interchange mode found in Ref. [23 |. 

Available data from the TORPEX device indicate that 
the magnitude of suprathermal ion dispersion is consis- 
tent with simulations at a single point in time during 
the ballistic phase [25|. Our estimates suggest that the 



transition from ballistic to interaction-phase values of 
should be observable for some experimental values of £. 
However, the experimental measurements of jr made so 
far are not sufficient to confirm our theoretical predic- 
tion, since a high resolution in the toroidal direction is re- 
quired, which will be possible to achieve with a toroidally 
moving source under construction. Measuring a change 
in jR into the asymmetric phase will be difficult because 
most of the ion beam tends to exit the plasma before this 
phase is well resolved. 

The interplay of fundamental phenomena such as gy- 
romotion and curvature drift, which determine the trans- 
port of suprathermal ions in the SMT, is also present in 
fusion confinement configurations. For example, orbit- 
averaging for trapped particles in a tokamak is analogous 
to drift-averaging for suprathermal ions in TORPEX. 
Also, radial constraint of transport due to zonal flows 
associated with ion-temperature gradient (ITG) turbu- 
lence [2^ or other velocity-shearing mechanisms may 
lead to similar subdiffusive tendencies as we find to be 
caused by vertical drift in the SMT. Our results cover 
10 kj_p < 1 to k±p - 10, and Kubo num- 

bers [23, nil in the range 0.2 < K. = v±Tc/Xc < 5. Here, 
Tc and Ac are the correlation time and correlation length 
of the turbulence, respectively. These ranges are relevant 
for present fusion devices with neutral beam injection Q 
and may be relevant for a future DEMO tokamak as well. 
While the SMT includes much of the fundamental physics 
for fusion plasmas, the general rules found here are rel- 
evant in other contexts as well, such as cosmic rays and 
solar flares 2^ 3^ . 
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